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Abstract 

Background: As the global population continues to expand, increasing yield in bread wheat is of critical 
importance as 20% of the world's food supply is sourced from this cereal. Several recent studies of the molecular 
basis of grain yield indicate that the cytokinins are a key factor in determining grain yield. In this study, cytokinin 
gene family members in bread wheat were isolated from four multigene families which regulate cytokinin synthesis 
and metabolism, the isopentenyl transferases {IPT), cytokinin oxidases (OCX), zeatin O-glucosyltransferases (ZOG), 
and (3-glucosidases {GLU). As bread wheat is hexaploid, each gene family is also likely to be represented on the A, 
B and D genomes. By using a novel strategy of qRT-PCR with locus-specific primers shared among the three 
homoeologues of each family member, detailed expression profiles are provided of family members of these 
multigene families expressed during leaf, spike and seed development. 

Results: The expression patterns of individual members of the IPT, CKX, ZOG, and GLU multigene families in wheat 
are shown to be tissue- and developmentally-specific. For instance, TalPT2 and TaCKXl were the most highly 
expressed family members during early seed development, with relative expression levels of up to 90- and 900-fold 
higher, respectively, than those in the lowest expressed samples. The expression of two cis-ZOG genes was sharply 
increased in older leaves, while an extremely high mRNA level of TaGLUl-1 was detected in young leaves. 

Conclusions: Key genes with tissue- and developmentally-specific expression have been identified which would be 
prime targets for genetic manipulation towards yield improvement in bread wheat breeding programmes, utilising 
TILLING and MAS strategies. 



Background 

Bread wheat (Triticum aestivum L.) is one of the worlds 
most important crop plants with a total production of 
over 600 million tons in 2007 [1]. This accounts for over 
20% of the worlds food supply. The increase in the 
global population is expected to increase the demand for 
wheat to over 800 million tons in 2030 and to more than 
900 million tons by 2050 [2]. Meeting this demand is 
predicted to depend largely on increasing yield through 
genetic improvement. However, in recent decades, inter- 
nationally coordinated public wheat breeding efforts have 
focused on increasing resistance to disease and abiotic 
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stress, while raising the genetic yield potential per se 
has received scant attention. Moreover, the challenges of 
increasing the genetic yield potential are considerable [3]. 

Grain size and grain number are significant compo- 
nents of yield in cereals, and both are major contributors 
to seed sink strength [4]. Yet, until recently, our under- 
standing of yield at the molecular level was limited and 
little progress has been made towards the identification 
of the genetic components that define yield in a quanti- 
tative manner [5]. Some of the most significant recent 
studies towards the understanding of the molecular basis 
of grain size and number include identification of a 
number of genes and quantitative trait loci (QTLs) con- 
tributing to seed number and size in Arabidopsis [6,7], 
rice [8,9] and maize [10], and the recognition that these 
genes and/or QTLs are cytokinin synthesis and metabol- 
ism genes. 
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Cytokinins (CKs) are hormones that regulate many 
developmental and physiological processes in plants 
[11], including reproductive development and grain yield 
[7,8,12,13]. In general, cytokinins may affect seed yield 
by enhancing sink strength and, consequently, increasing 
seed size, by up-regulating cell cycle related genes, 
through the action of sugar signalling and by enhancing 
phloem unloading and sugar import into the endosperm 
through enzymatic activity of cell wall invertase [12]. 
Although the beneficial effect of exogenously applied 
cytokinins on grain yield has been documented for 
several decades, including for wheat, maize, rice and 
barley [14-17], the application of cytokinins to field 
crops is limited due to the substantial cost of cytoki- 
nins. Given that changes in the levels of endogenous 
cytokinins will alter the regulation of many physio- 
logical processes resulting in the disruption of normal 
growth patterns [18], and rapid and dramatic changes 
in cytokinin concentration have been observed during 
panicle and early grain development in a number of 
crop species including wheat [19-21] and rice [8,22], 
yield improvement through genetic manipulation of 
cytokinin synthesis and metabolism genes to artificially 
manipulate the endogenous cytokinin levels could be of 
considerable importance. 

The last decade has seen the key steps in cytokinin 
biosynthesis, metabolism, and signal transduction eluci- 
dated at the molecular level [23]. The bioactive cytokinin 
concentration is regulated by several multigene families, 
including isopentenyl transferases (IPT) for cytokinin 
biosynthesis, cytokinin oxidases (CKX) for degradation, 
zeatin O-glucosyltransferases (ZOG) for reversible inacti- 
vation, and p-glucosidases (GLU) for reactivation. In 
Arabidopsis, nine IPT genes, seven CKX genes, three 
ZOG genes and GLU genes have been identified and 
functionally verified [24-28]. Their orthologues have also 
been annotated and/or functionally verified in a number 
of important crop species. In rice, for instance, eight 
IPTs and 11 CKXs have been isolated [29-31]. Similarly, 
11 IPTs, 13 CKXs, and several ZOGs and GLUs have 
been isolated from maize [32-34]. 

Because of the multigene family nature of these cytoki- 
nin metabolism genes, each member can be functionally 
differentiated based on their spatial and temporal 
expression patterns [31]. Indeed, the temporal and spatial 
expression, and functional analysis, of specific family 
members has been confirmed in a number of species. 
For example, in Zea mays, the temporal and spatial 
expression, and functional analysis, of specific gene family 
members, confirmed that maize seed biosynthesises its 
own endogenous cytokinin, and that ZmIPT2 and 
ZmCKXl are the gene family members controlling cytoki- 
nin homeostasis in the developing maize grain [32,35]. 
In rice, Ashikari et al [8] showed that a QTL associated 



with increased grain productivity was in fact a mutated 
cytokinin oxidase gene, OsCKX2. Most recently, RNAi 
was used to silence HvCKXl in barley with a resultant 
increase in both grain number and grain weight [13]. 

Due to the pleiotropic nature of cytokinins and the 
variation in spatial expression of gene family members, 
manipulation of a particular gene may cause side effects, 
affecting non target organs or developmental stages. 
Therefore, detailed expression profiles of all available 
members of the major gene families involved in cytoki- 
nin synthesis and metabolism (i.e. the IPTs, CKXs, ZOGs 
and GLUs), should be investigated simultaneously in 
order to determine the key genes with critical roles in the 
target organ or developmental stage. While there are 
examples in the literature showing that there is variation 
in the spatial and temporal expression patterns of individ- 
ual cytokinin gene family members [32,36-41], simultan- 
eous monitoring of the four major cytokinin metabolism 
gene families has only recently been reported in maize 
[33], but without detailed information during leaf, floral 
and seed ontology. 

Apart from our preliminary report [41], detailed infor- 
mation regarding the expression profiles and functional 
territories of these multigene families is lacking for bread 
wheat, possibly due to the added complexity of its hexa- 
ploid nature, such that each gene is represented in tripli- 
cate homoeologues on the three genomes [42]. Further, 
each triplicate homoeologue within the same locus can 
be differentially regulated by genetic and epigenetic 
mechanisms in a tissue and developmentally specific 
manner [43-45]. In general, uniformly expressed homo- 
eologues only account for 20-40% of the loci in bread 
wheat, with 20-30% of loci presenting at least one 
silenced homoeologue and, for the rest of the loci, the 
three homoeologues are differentially expressed across 
organs and developmental stages [43-46]. 

Given the above, we hypothesised that cytokinin 
homeostasis within an organ is co-ordinately regulated 
by particular homoeologues of different gene family 
members to allow the precise control of organ develop- 
ment, and that optimum seed yield could be obtained 
by manipulating the endogenous cytokinin concentra- 
tion at specific stages of reproductive growth and seed 
development through genetic manipulation utilising, for 
example, over-expression, gene silencing, TILLING and 
MAS techniques. 

However, it would be very complicated and essentially 
impractical to determine the expression profile of each 
homoeologue of each gene family member simultan- 
eously, particularly when multiple tissue types and time 
points are involved. Due to this complexity only three 
TaCKX genes [47-49] and one set of GLU homoeologous 
genes [50] have been isolated and investigated so far, 
with little information regarding their expression profiles. 
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No study with regard to /PTand ZOG genes has yet been 
reported for wheat. 

In this study, we report the isolation and simultaneous 
expression analysis of TalPT, TaCKX, TaZOG and 
TaGLU gene family members during flag leaf, spike and 
seed development by using a novel strategy in which 
locus-specific PCR primers shared by the three homoeo- 
logues were used, and reveal the co-ordinated expression 
of several members of different multigene families. 

Methods 

Plant material 

Winter bread wheat, variety Equinox, was grown under 
prevailing climatic conditions in a nursery plot at the 
New Zealand Institute for Plant & Food Research, Lin- 
coln, Canterbury, New Zealand. After microscopic deter- 
mination, representative developmental stages of spikes, 
carpels, seeds and the basal one-third of flag leaves were 
harvested, from spikelet initiation to 28 days after anthe- 
sis (daa) and immediately frozen in liquid nitrogen and 
stored at -80 °C until used. The chlorophyll content of 
flag leaf samples was determined as described in O'Keefe 
et al. [40]. The total chlorophyll content was calculated 
using the equations of Porra [51]: 

Chi a + b = 17.67Abs 647 + 7.12Abs 664 . 

RNA extraction and cDNA synthesis 

Total RNA was extracted using TRIzol Reagent (Invi- 
trogen, Carlsbad, CA, USA) following the manufacturer s 
protocol. Two independent tissue samples of each devel- 
opmental stage were used as biological replicates. The 
concentration and purity of the total RNA extracted 
was determined using a Nanodrop spectrophotometer 
(Nanodrop Technologies Inc, Wilmington, DE, USA). 
The integrity of the total RNA was determined by elec- 
trophoresis on 1% (w/v) agarose gels. For cDNA synthe- 
sis, 1 \xg total RNA, 50 U Expand Reverse Transcriptase 
(Roche, Mannheim, Germany), 50 pmol oligo (dT) pri- 
mers and 100 pmol random hexamer (pdN6) primers 
were used in a 20 \A reverse transcription reaction. The 
final reaction mix was incubated at 42 °C for 2 h. 

Gene isolation and sequence analysis 

To identify putative IPT, CKX, ZOG and GLU ortholo- 
gue sequences in wheat, all the annotated family mem- 
bers of these multigene families from Arabidopsis, rice, 
maize and barley were used as query sequences to 
BLAST search the GenBank database (http://blast.ncbi. 
nlm.nih.gov). Subsequent verification utilised the recently 
released draft wheat genome (http://www.cerealsdb.uk. 
net/search_reads.htm). Isolation of these putative ortho- 
logues in wheat was carried out through direct sequen- 
cing of RT-PCR products using cDNA template derived 



from a mixed tissue sample containing several stages of 
developing leaves, spikes, carpels and seeds. PCR primers 
specific to each family member were designed based on 
the obtained sequence information. For family member (s) 
without sequence information in the available databases, 
degenerate primers were designed within the conserved 
regions of the same family member in barley, rice, maize 
and other monocot species, using Primer Premier 5.0 
(Additional file 1). 

PCR amplification was conducted in a 20 \A reaction 
mix containing 2 mM MgCl 2 , 10 pmol of each forward 
and reverse primer, 2 \A of 10-fold diluted cDNA and 
1 U FastStart Taq polymerase (Roche), using cycling 
conditions as described in Song et al. [52]. The PCR 
products were separated on a 1% agarose gel. All bands 
of approximately the expected size were purified and 
sequenced on an ABS 3730 sequencer (Applied Bio- 
science) using the same primers as for the PCR. The raw 
sequences from at least two forward and two reverse 
directions were aligned using ClustalX software to cor- 
rect the sequencing errors and to generate a consensus 
sequence for each gene of interest (GOI). The sequences 
were verified by BLAST searching the GenBank data- 
base. Full length or 3' end cDNA sequences of selected 
genes were generated using gene specific primers and 
the GeneRacer™ Kit or 3' Race Kit (Invitrogen) follow- 
ing the manufacturer s protocol. For phylogenetic ana- 
lysis, orthologues from Arabidopsis, maize, rice and 
other closely related monocot species for each GOI were 
used to construct both a Neighbor-Joining (NJ) phylo- 
gentic tree using ClustalX software [53] and a Maximum 
Likelihood (ML) tree using PhyML3.0 software [54,55], 
with 1000 bootstrap replicates. Each tree was rooted 
with an out group orthologue. The nucleotide sequences 
reported in this paper have been submitted to the Gen- 
Bank under the accession numbers listed in Table 1. 

Quantitative expression of cytokinin synthesis 
and metabolism genes using qRT-PCR 

Due to the hexaploid nature of bread wheat, in which 
each gene is present in triplicate sets located on homolo- 
gous chromosomes of genomes A, B and D [42], gene 
expression studies involving large numbers of genes, 
particularly those in multigene families, is very challen- 
ging. Since, in general, the three homoeologues of the 
same locus share high sequence similarity [56], we 
employed a novel strategy to minimize the work load 
and to handle, simultaneously, as many family members 
as possible. PCR primers specific to the family member 
but which can be shared by the three homoeologues of 
the same gene were designed and used for quantitative 
expression studies using qRT-PCR (Additional file 2). 
The temporal expression of each putative triplicate gene 
set and of selected reference genes was quantified using 
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Table 1 Sequence similarity of the genes of interest (GOIs) isolated in the present study with their othologues 



GOI in this study 



Orthologues in closely related species 



Similarity (%) 



Gene name 


Accession No. 


Gene name 


Accession No.* 


Species 


Protein 


Nucleotide 


TolPT2 


JN 128577 


ZmlPT2 


ABY78882 


Zea mays 


76.8 


79.9 


TalPT3 


JN 128578 


OslPT8 


AB239805 


Oryza sativa 


84.9 


85.5 


TalPT5 


JN 128579 


ZmlPT5 


NP001 1 21 1 94 


Zea mays 


78.1 


82.0 


T0IPT6 


JN 128580 


ZmlPT6 


ABY78885 


Zea mays 


80.4 


83.9 


TalPT7 


JN128581 


ZmlPT7 


ABY78886 


Zea mays 


79.8 


82.0 


T0IPT8 


JN 128582 


ZmlPT8 


ABY78887 


Zea mays 


80.4 


82.8 


ToCKXl 


JN 128583 


ZmCKXl 


NP001 105591 


Zea mays 


75.1 


84.4 


ToCKX2 


JN 128584 


ZmCKX5 


NP001 185958 


Zea mays 


68.0 


78.0 


TaCKX3 


JN 128585 


ZmCKXIO 


NP001 146838 


Zea mays 


79.2 


84.9 


TaCKX4 


JN 128586 


OsCKX4 


NP001 045353 


Oryza sativa 


88.9 


76.1 


T0CKX6 


JN 128587 


OsCKX3 


N PO0 1064886.1 


Oryza sativa 


88.3 


76.1 


ToCKX7 


JN 128588 


OsCKX6 


Q6YW51 


Oryza sativa 


72.7 


75.6 


TaCKX8 


JN 128589 


ZmCKX8 


NP001 185809 


Zea mays 


73.5 


76.3 


TaCKX9 


JN 128590 


OsCKXIO 


Q5Z620 


Oryza sativa 


72.9 


71.9 


TaCKXIO 


JN128591 


OsCKX9 


Q75K78 


Oryza sativa 


86.9 


72.4 


TaCKXl 1 


JN 128592 


OsCKX8 


A2XVN3 


Oryza sativa 


81.5 


74.6 


TacZOGI 


JN 128593 


ZmcZOGl 


NP001 105017 


Zea mays 


75.5 


82.0 


TacZ0G2-1 


JN 128594 


ZmcZOG2 


Q8RXA5 


Zea mays 


75.4 


82.6 


TocZOG2-2 


JN 128595 


ZmcZOG2 


Q8RXA5 


Zea mays 


75.2 


80.3 


ToZOGl 


JN 128596 


ZmZOG3 


NO001 148195 


Zea mays 


74.9 


82.2 


TaZ0G2 


JN 128597 


ZmcZOG2 


NP001 148091 


Zea mays 


64.1 


75.6 


TaZ0G3 


JN 128598 


ZmZOG3 


NO001 148195 


Zea mays 


72.3 


80.7 


ToGLUl-1 


JN 128599 


ZmGLU2 


ACG24271 


Zea mays 


67.9 


73.3 


TaGLUl-2 


JN 128600 


ZmGLU2 


ACG24271 


Zea mays 


67.3 


73.1 


ToGLUl-3 


JN 128601 


ZmGLU2 


ACG24271 


Zea mays 


68.0 


74.1 


ToGLU2 


JN 128602 


ZmGLUl 


NM001 1 1 1984 


Zea mays 


95.2 


94.9 


ToGLU3 


JN 128603 


OsGLU25 


Q0DA21 


Oryza sativa 


85.8 


NV 


ToGLU4 


JN 128604 


HvGLU 


ACF07988 


Hordeum vulgare 


94.8 


93.6 


TaRRl 


JN 128605 


OsRR1 


CAI79405 


Oryza sativa 


73.7 


79.7 


TaRR4 


JN 128606 


OsRR4 


N POO 1045420 


Oryza sativa 


90.0 


84.9 


TaRR5 


JN 128607 


OsRR4 


NP001 053348 


Oryza sativa 


95.0 


89.4 


ToRR9 


JN 128608 


OsRR9 


NP001 065722 


Oryza sativa 


82.4 


82.6 



* If accessions numbers are different for protein and nucleotide sequences of the same gene, those for protein sequences are listed. NV, corresponding sequence 
not available. 



an Mx3000P R real-time PCR instrument (Stratagene) 
and a home-made SYBR Green master mix or the KAPA 
SYBR® FAST qPCR Kits (Kapa Biosystems, Boston, 
USA). PCR amplification efficiency of each GOI and 
reference genes were determined as described in Song 
et al. [57]. At least two technical replicates for each of 
two biological replicates were carried out for samples 
from each developmental stage 

Because the changes in endogenous cytokinins immedi- 
ately post anthesis are well recorded [19-22], we thought 
to use type A response regulators as a surrogate indicator 



of endogenous cytokinin levels. These genes have recently 
been used as markers to track the in planta bioactive 
cytokinin concentration because their expression levels 
are reported to be proportionally correlated with the 
endogenous cytokinin content of the tissue [58,59]. 

Multiple reference genes, GAPDH, fi-actin, 18 S rRNA, 
and protein phosphatase gene (PP2A), were used as 
internal controls. For each reference gene, a correction 
factor (CF) for each cDNA sample was calculated using 
the Ct number of this cDNA sample divided by the aver- 
age Ct number of all cDNA samples included in the 
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same experiment. The values of three to four technical 
replicates were averaged to form the CF for each bio- 
logical replicate of each reference gene. The final CF value 
for each biological replicate was created by averaging the 
CF values of the four reference genes. For each cDNA 
sample, the Ct number of each target gene was corrected 
before statistical analysis of its expression level. 

Data analysis 

Experimental data were analysed as a randomised block 
experiment by analysis of variance as described in Song 
et al [52], where PCR-run was the block factor, the stage 
of development was the treatment factor, and relative 
transcript abundance was the response variable. 

Results 

Isolation of cytokinin synthesis and metabolism genes 

After using all of the annotated family members of IPT, 
CKX, ZOG and GLU in Arabidopsis, rice, maize and bar- 
ley as query sequences to BLAST search the publically 
available genomic databases, we identified a large num- 
ber of putative IPT, CKX, ZOG and GLU orthologue 
sequences of these multigene families in wheat. Most of 
the sequences were ESTs from GenBank and genomic 
DNA fragments from the recently released draft wheat 
genome (http://www.cerealsdb.uk.net/search_reads.htm). 
After manual assembly of the genomic fragments, all 
putative sequences were multiple aligned and phylogen- 
etically analysed to determine the members of each 
gene family. 

For the majority of members of these multigene fam- 
ilies, three groups of sequences with very high similarity 
were detected. They were considered to be the three 
homoeologues of a putative gene. At least one set of spe- 
cific PCR primers was designed for RT-PCR isolation of 
this gene. In total, we have isolated at least one homoeo- 
logous sequence from six IPTs, ten CKXs, five ZOGs and 
four GLUs (Table 1). Phylogenetic analysis (both ML 
and NJ trees) grouped each gene family into several dis- 
tinct clusters consisting of gene sequences from both 
dicot and monocot species, with high bootstrap confi- 
dences between monocot and dicot sub-clusters. The 
bootstrap confidence levels among clusters were also 
high for IPT, ZOG and GLU genes while those for CKX 
genes were low (Figures 1, 2, 3, 4. Additional files 3, 4, 5, 
6). Multiple sequence alignment results showed that 
CKX genes from different groups shared much greater 
sequence similarities than other gene families (Add- 
itional files 7, 8, 9, 10). Phylogenetic analysis also showed 
that these newly isolated sequences represented most of 
the family members in this species as predicted accord- 
ing to other closely related monocot species for which 
whole genome sequences were available, such as rice 
and maize (Figures 1, 2, 3, 4. Additional files 3, 4, 5, 6). 
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Figure 1 The maximum likelihood (ML) phylogenetic tree for 
IPT proteins in Arabidopsis thaliana (AtlPT), Oryza sativa 
(OslPT), Triticum aestivum (TalPT), and Zea mays (ZmlPT). The 

tree was rooted using IPT protein from Rhodococcus foscions (RflPT). 
Node values are percentages of bootstraps generated with 1000 
bootstrap replicates. 



For the IPT gene family, TaIPT2, 3, 5, 6, 7 and 8 were 
similar, respectively, to ZmIPT2, 3, 5, 6, 7 and 8 in maize 
and were designated TaIPT2, 3, 5, 6, 7 and 8, respect- 
ively (Table 1, Figure 1, Additional files 3 and 7). TaIPT2 
is a 1012 bp fragment, with deduced protein sequence of 
284 aa from the start point of the coding region and 
160 bp of the 5' UTR sequence. TalPT 5 and 6 frag- 
ments are 1561 and 1552 bp, respectively, with deduced 
protein sequences of 470 and 433 aa, respectively, from 
the start point of the coding region, and 151 and 253 bp 
of the 5' UTR sequences, respectively. TaIPT3, 7 and 
8 are 694, 753 and 950 bp fragments within the coding 
region of the genes. 
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Figure 2 The maximum likelihood (ML) phylogenetic tree for 
CKX proteins in Arabidopsis thaliana (AtCKX), Oryza sativa 
(OsCKX), Triticum aestivum (TaCKX), and Zea mays (ZmCKX). 

The tree was rooted using CKX protein from Rhodococcus fascians 
(RfCKX). Node values are percentages of bootstraps generated with 
1000 bootstrap replicates. 



For the CKX gene family, three sets of TaCKXs were 
almost identical to the recently annotated wheat 
TaCKXI, 2, and 3 [47-49] and were named accordingly. 
These were putative orthologues of ZmCKXI, 5 and 10, 
respectively, in maize. Other putative TaCKXs grouped 
together with ZmCKX4, 6-11, and were designated 
TaCKX4, 6-11, respectively (Table 1, Figure 2, Additional 
files 4 and 8). TaCKXI, 2, 3, 6, 7 and 8 are cDNA frag- 
ments of 1500-2000 bp covering the whole coding 
region and contain fragments of both the 3 ' and 5 ' UTR 



sequences. TaCKX4 and 9 fragments are 989 and 
788 bp, respectively, with deduced protein sequences 
of 199 and 252 aa, respectively, near the 3' end of the 
coding region, and 392 and 32 bp of the 3' UTR 
sequences, respectively. TaCKXIO and 11 are 541 and 
735 bp fragments within the coding region of the genes. 

Three putative ZOG genes shared high sequence simi- 
larity with ZmZOGl, ZmZOG2 and ZmZOG3, and were 
named TaZOGl, TaZOG2 and TaZOG3, respectively. 
Two other putative ZOG genes were orthologues of cis- 
ZOG1 and cisZOG2 in maize, and designated Tacis- 
ZOG1 and TacisZOG2 (Table 1, Figure 3, Additional 
files 5 and 9). TaZOGl and 3, and TacisZOGl and 2-1 
are cDNA fragments of around 1400 bp covering the 
whole coding region and contain short fragments of 
both the 3' and 5' UTR sequences. TaZOG2 and Tacis- 
ZOG2-2 are 1100 and 963 bp, respectively, with deduced 
protein sequences of 366 and 322 aa, respectively, within 
the coding region of the genes. 

In addition to the one triplicate set of GLU genes, 
TaGLUla, lb and lc, investigated by Sue et al [50], we 
have isolated three other GLU genes, TaGLU2, 3 and 4, 
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Figure 3 The maximum likelihood (ML) phylogenetic tree 
for ZOG proteins in Arabidopsis thaliana (AtZOG, UGT) 7 
Triticum aestivum (TaZOG, TacZOG), Zea mays (ZmZOG, 
ZmcZOG) and Phaseolus lunatus (ZOG1). The tree was rooted 
using cytokinin-N-glucosyltransferase 1 protein from Zea mays 
(ZmZNGI). Node values are percentages of bootstraps generated 
with 1000 bootstrap replicates. 



Song et al. BMC Plant Biology 201 2, 12:78 
http://www.biomedcentral.com/1471-2229/12/78 



Page 7 of 16 



100 



100 



48.2 



BnBGL 

AtBGL19 

TaGlu3 

OsGki25 



98.9 



97.6 



TaGlu4 
HvACF07998 



98.3 



95,1 



58,1 



OSB7F7K7 

ZmGlul 

TaGlu2 



69.5 



70.5 



ZmGlu2 

TaGlulb 

TaGiula 



TaGlulc 



OsBGal3 

Figure 4 The maximum likelihood (ML) phylogenetic tree 
for representative GLU proteins. At, Arabidopsis thaliana; Os, 
Oryzo sotiva; Ta, Triticum oestivum; Zm, Zeo mays. The tree was 
rooted using (3-galactosidase 3 protein from Oryzo sativa (OsBGaB). 
Node values are percentages of bootstraps generated with 1000 
bootstrap replicates. 



which shared high sequence similarity to ZmGLU2 in 
maize, and OsGLU25 and B7F7K7 genes in rice, respect- 
ively (Table 1, Figure 4, Additional files 6 and 10). 
TaGLUl-1, 1-2, 1-3 (the equivalents of TaGLUla, lb 
and lc [50]), TaGLU2 and TaGLU4 are cDNA fragments 
of 946-1648 bp in length near the 3' end of the coding 
region, with deduced protein sequences of 277-501 aa, 
and 110-361 bp of the 3' UTR sequences. TaGLU3 is a 
1080 bp fragment at the 5 ' end of the coding region. 



Six type-A response regulators, TaRRl, 3-6 and 9 were 
isolated. They were putative orthologues of OsRRl, 3-6 
and 9, respectively (Table 1). 

Chlorophyll content during flag leaf development 
and senescence 

The total chlorophyll content of the basal one -third of 
the flag leaves at different developmental stages was 
determined (Figure 5). The chlorophyll content in the 
flag leaf increased rapidly as the young leaf expanded to 
its full size, remained high (over 4 mg/g FW) over the 
next four weeks, from fully expanded to 7 daa, and then 
started to decrease. At 14 daa, when the end one-third 
of the leaf blade showed signs of senescence and the rest 
of the leaf remained green to the eye, the chlorophyll 
content in the basal one-third of the leaf had dropped to 
around 3.6 mg/g FW. At 21 daa, when the end two-fifths 
of the leaf were yellow in colour, about half of the 
chlorophyll in the basal one-third of the leaf was 
degraded, although visually the colour of this part of the 
leaf still remained green. At 28 daa, half of the leaf blade 
was yellow and the end one third had started to dry out. 
At that time, over 80% of the chlorophyll was degraded 
in the basal one-third of the leaf blade, where significant 
signs of senescence were now observable. 

Quantitative expression of cytokinin synthesis 
and metabolism genes 

In order to determine the key genes with likely/critical 
roles during flag leaf, spike, carpel and early seed devel- 
opment, we first did a preliminary screen of all available 
TalPT, TaCKX, TaZOG and TaGLU members utilizing 
cDNAs from root tips, elongating zone of internodes, 
leaves and reproductive tissues of mixed developmental 
stages. Gene specific PCR primers targeted to the three 
homoeologues of the same locus were used for qRT- 
PCR. Only those family members which expressed 
exclusively or mainly in leaves and reproductive tissues 




2 cm 15 cm Expd 7dba 1 daa 7 daa 14 daa 21 daa 28 daa 



Figure 5 Total chlorophyll content in flag leaf samples used for mRNA quantification of cytokinin regulatory genes in wheat. Flag leaf 
sample order: 2 and 15 cm in length, expanded, days before and after anthesis. 
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were then selected for further detailed expression studies 
across an extended series of spike, carpel, seed and flag 
leaf developmental stages. SEM and light microscopy 
were used to record the precise spike, carpel and seed 
developmental stages of the samples (Figure 6). 

Using GAPDH, fi-actin, 18 S rRNA, and PP2A as refer- 
ence genes, marked changes in the expression of mem- 
bers of all four gene families and of several type-A 
response regulators were observed during reproductive 
and flag leaf development (Figures 7 and 8). 

During reproductive development, TaIPT2 was highly 
expressed during early seed development (Figure 7A). Its 
relative expression level was low during early develop- 
ment of the spike and carpels but increased dramatically 
1 to 4 daa, being 60- to 90-fold higher than the baseline 
expression, and 5 to 7-fold higher than that at 1 day 



before anthesis (dba). The mRNA abundance of TaIPT2 
dropped sharply 7 daa and returned to near the baseline 
at 21 daa. Both TalPTS and TaIPT8 showed three peaks 
of expression during reproductive development. One 
peak was shared by both genes and occurred around 
2 daa, but with much less mRNA abundance than that 
of TaIPT2. While the other two expression peaks of 
TalPTS appeared in 10 cm spikes and mature ovules at 
1 dba, respectively, those for TaIPT8 appeared in young 
spikes of less than 2 cm in length and mature carpels at 
7-14 dba. The expression of TalPT genes was, in gen- 
eral, lower during leaf development than reproductive 
development. While an expression peak appeared in the 
young leaf for TaIPT2 and in the expanded leaf for TalPTS, 
the two genes shared another expression peak at 7 daa. In 
comparison with reproductive development, the expression 
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Figure 6 Developmental stages of spike, carpel and seed samples used for mRNA quantification of cytokinin regulatory genes in 
wheat observed using SEM (A) and light microscopy (B-L). A, spike 0.5 cm in length; B, spike 2 cm in length; C, spike 10 cm in length; 
D, carpel 14 days before anthesis (dba); E, carpel 7 dba; F, carpel 1 dba; G, 1 day after anthesis (daa); H, 2 daa; I, 4 daa; J, 7 daa; K, 14 daa; L, 
21 daa. a, anther; c, carpel; e, endosperm; fl-f5, florets 1-5 in a spikelet; fm, floral meristem; slm, spikelet meristem. Scale bars: A = 50 pm, 
B-E = 500 pm, F-K= 1 mm. 
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Figure 7 (See legend on next page.) 
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(See figure on previous page.) 

Figure 7 Quantitative expression profiles of selected putative cytokinin regulatory genes during carpel and seed development in 
bread wheat. A, TolPT genes; B, TaCKX genes; C, ToZOG genes; D, TaGLU genes; E, ToRR genes. Data values are means of relative mRNA levels in 
fold changes detected using qRT-PCR, using at least two technical replicates for each of the two biological replicates. Error bars represent the SD 
calculated for the combined technical and biological replicates. GAPDH, (5-actin, 18 S rRNA, and protein phosphatase gene {PP2A), were used as 
internal controls. Before quantification of the expression level of each of the target genes, the Ct numbers for each target gene were corrected 
by using the average correction factor determined for each of the four reference genes. 



of TaIPT8 was very low during all leaf development stages. 
The TaIPT3 and TaIPT7 genes expressed at very low levels 
in all samples (Figure 8A). 

For the TaCKX gene family, TaCKX 1 was the most 
highly expressed family member, with expression from 
before anthesis through to 14 daa and relative expression 
levels of 220 to 980-fold higher than that in the lowest 
expressed sample (Figure 7B). The expression pattern of 
TaCKX2 was similar to that of TaCKX 1, apart from the 
much lower expression level and the postponed expres- 
sion window. While TaCKX6 expressed during all stages 
of reproductive development, only one expression peak 
was detected at 7dba. In leaf samples, TaCKXl was the 
only family member with high expression, particularly in 
mature and senescing leaf tissue starting around 14 daa, 
when the chlorophyll content in the basal one- third just 
started to decrease. TaCKX2 expressed constantly across 
all leaf samples but at a relatively low level (Figure 8B). 
Expression of the other CKX genes, i.e. TaCKX3, 4 and 
7, remained very low during both reproductive and leaf 
development, although slightly higher mRNA levels were 
detected in 10 cm spikes and young leaves for TaCKX4, 
and TaCKX3, respectively (Figures 7B, 8B). 

In comparison with the senescing flag leaf, zeatin O- 
glucosyltransferase (ZOG) genes expressed at relatively 
low levels during reproductive development (Figure 7C). 
However, expression of TacisZOGl showed a peak of 
activity during spike development, as well as expression 
from 1 dba and throughout early grain development. 
Additionally, TacisZOG2-2 expressed from 1 to 14 daa. 
Both TacisZOGl and TacisZOG2-l expressed during leaf 
expansion and in the mature leaf until anthesis. In older 
leaves, while other ZOG genes expressed at very low 
levels, there was a sharp increase of mRNA for both 
TacisZOGl and TacisZOG2-l, The expression of Tacis- 
ZOGl was over 900-fold higher than that in the lowest 
expressed sample tested and was particularly high from 
21 daa, when the chlorophyll content in the basal one- 
third of the leaf started to decrease. Expression of one 
homoeologue of TacisZOG2, TacisZOG2-l mimicked the 
trend of TacisZOGl in leaves, whereas the TacisZOG2-2 
homoeologue showed virtually no expression (Figure 8C). 
In contrast, in the developing ovule and young seed, 
while the TacisZOG2-2 homoeologue was showing a 
strong pattern of expression, the TacisZOG2-l showed 
only low, fluctuating expression. 



For the /?-glucosidase (GLU) gene family, an extremely 
high expression of TaGLU 1-1, with a relative expression 
level over 1400-fold, was detected in young expanding 
leaves 2 cm in length (Figure 7D). Multiple expression 
peaks of the same gene, with relative expression over 
200- to 800-fold, were also detected in 0.5 and 10 cm 
spikes, and in leaves at anthesis. In contrast, another al- 
lele from the same homoeologue gene, TaGLU 1-2, 
expressed consistently during spike, seed and leaf devel- 
opment, with the exception of in leaves at anthesis, 
during which the relative expression level was 300-fold 
higher than the baseline and similar to TaGLUl-1. The 
highest expression of TaGLU4 was detected at 7 dba, 
presenting an over 500-fold up regulation. TaGLU4 also 
expressed at a relatively high level at most other stages 
of reproductive development, but at a consistently lower 
level in most leaf samples. Expression of other GLU fam- 
ily members, such as TaGLU2 and TaGLU3, was low 
across all tested samples, without any obvious expression 
patterns during reproductive and leaf development 
(Figure 8D). 

With the exception at the time of anthesis, elevated 
expression of several response regulators (RR) was 
detected in fully expanded and mature leaves (Figure 8E). 
This expression was coincident with that of TalPT gene 
expression in the leaves. Although several RR members 
did express at relatively high levels in the early stages of 
spike development (Figure 7E), there was only a rela- 
tively low level of RR expression during reproductive 
development when TalPT gene family members were 
expressing. However, the small peak of TaRR9 is coinci- 
dent with the peak of IPT expression at 2 daa. 

Discussion 

Phylogenetic analysis showed that TaCKXl shares high 
sequence similarity with orthologues from other monocot 
species, including HvCKXl in barley, ZmCKXl in maize 
and OsCKXl in rice (Figure 2), suggesting that compara- 
tive phylogenetic analysis based on sequence similarity is 
a powerful method to identify homologous genes and to 
predict their physiological functions in target species. 
Identification of orthologues across species is essential to 
predict gene function in newly sequenced genomes such 
as wheat for which a draft genome has recently been 
released (http://www.cerealsdb.uk.net/search_reads.htm). 
Indeed, when Gu et al [34] combined a phylogenetic 
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Figure 8 Quantitative expression profiles of selected putative cytokinin regulatory genes during flag leaf development in bread wheat. 

See Figure 7 for legend. 
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analysis with comparative expression of Arabidopsis, pop- 
lar, maize and rice OCXs in various tissues, they observed 
a general association of sequence similarity with expres- 
sion patterns. Phylogenetic analysis has been widely used 
to classify gene families and predict their functional ortho- 
logues, and has been successfully used to generate func- 
tional categories within large gene families such as the 
WRKY and bHLH transcription factor families [60,61]. 
However, phylogenetic analysis does not allow the differen- 
tial expression of homoeologous genes to be determined. 

In this study, we isolated most members of the major 
multi-gene families, TalPTs, TaCKXs, TaZOGs, and 
TaGLUs involved in cytokinin synthesis and metabolism. 
Potential functional homoeologues from the three gen- 
omes, A, B and D, were pre-screened based on their 
expression in selected tissue types and in developmental 
stages closely related to grain yield determination. By 
using a novel strategy with specific PCR primers shared 
among the three homoeologues of each family member, 
we have been able to reveal detailed expression profiles 
of family members of the four multigene families during 
leaf, spike and seed development using qRT-PCR. Such 
information is required before determining which are 
the key genes and homoeologues contributing to cytoki- 
nin homeostasis within specific organs and developmen- 
tal stages. These key genes then become the target for 
further investigation and genetic manipulation toward 
yield improvement in a breeding programme. 

In plants, the rate limiting enzyme in cytokinin biosyn- 
thesis is considered to be IPT, which functions by 
attaching the isopentenyl side chain to the N 6 moiety of 
ADP or ATP. At least three of the six IP T genes tested 
in this study, TaIPT2, TalPTS and TaIPT8 showed spe- 
cific expressed patterns during spike, carpel, seed and 
leaf development. All three genes were expressed at the 
highest levels at early seed development immediately 
after anthesis (Figure 7A), a developmental window coin- 
cident with the cell division period, during which the cell 
number of the endosperm is determined, and during 
which endogenous cytokinin levels peak [e.g. 19-21]. 
Expression of IPT family members was also detected in 
young spikes and developing carpels, suggesting that 
these genes should be further investigated as players 
in determining the number as well as the size of the 
reproductive organs. Interestingly, these TalPT genes 
expressed at relatively low levels during the expansion 
and the fully expanded phase of the leaf. Together with 
similar results recently reported in maize [32,33], this 
suggests that it would be possible to manipulate IPT 
activity during reproductive development without signifi- 
cantly affecting normal leaf functions. An initial target 
for enhanced expression could be TaIPT8, which does 
not express highly in leaves, thereby avoiding creating 
competition between sink and source. 



The CKX gene family members have been shown to 
express in different plant tissues and to play essential 
roles in controlling cytokinin levels during plant growth 
and development [62]. Information about the expression 
and physiological function of TaCKXs was only recently 
available and only for TaCKXl, 2 and 3 [47-49]. While 
investigating the expression of ten sets of TaCKX genes, 
we showed that at least three genes are specifically 
or preferentially expressed during seed development 
(Figure 7B). The strong and specific expression of the 
TaCKXl gene during early seed development provides a 
useful candidate target. By perturbing cytokinin homeo- 
stasis during this key yield determining stage, i.e. by 
down regulating the expression of any one of CKX1 } 2 
or 6 during seed development and disturbing the normal 
mechanism of cytokinin homeostasis, the resulting ele- 
vated cytokinin level may lead to enhanced seed number 
and size [8]. Indeed, Zalewski et al [13] were able to 
increase grain yield of barley by down regulating the 
TaCKXl orthologue, HvCKXl, using RNAi. 

O-glucosylation is a major step in the metabolism of 
cytokinins. Since O-glucosylation is reversible and O-glu- 
cosides are resistant to cleavage of the A^-side chain by 
CKX, this conversion forms part of the homeostatic con- 
trol of active cytokinin levels through temporary storage of 
cytokinins [11]. Interestingly, TacisZOGl showed elevated 
expression pre-anthesis (Figure 3), at a time when en- 
dogenous O-glucosides start to accumulate in the develop- 
ing ovule [19]. Extremely high level of expression was 
detected for two cis-ZOG genes, TacisZOGl and Tacis- 
ZOG2-1. The dramatic increase of TacisZOG expression 
during leaf senescence in the present study agrees with 
early observations of high levels of O-glucosides in sen- 
escing leaves [63] and cytokinin inactivation by O-gluco- 
sylation promoted by senescence-related processes in 
petunia [64]. In contrast, Vyroubalova et al [33] only 
detected low expression of ZmcisZOG in old leaves of 
maize. This difference might be due to differences in 
cytokinin metabolism during senescence in these species 
[63]. 

The /?-D-glucosidases are members of the GH1 family, 
which is implicated in the hydrolysis of a number of 
plant secondary metabolites, including the regulation 
of the biological activity of cytokinins [65]. Until now, 
Zm-p60 1 ZmGLUl in maize, TaGLUl in wheat and Bgl4 
in Brassica napus are the only enzymes for which the 
ability to release cytokinins from O-glucoside conjuga- 
tion has been demonstrated [50,65-67]. In wheat, only 
one homoeologous set of GLU genes, TaGLUla, b and c, 
has been investigated [50], but with no information shown 
on in planta expression during reproductive and leaf 
development. Our qRT-PCR data showed high expression 
of TaGLUl-1, the equivalent of TaGLUla in Sue et al. 
[50], in 0.5 cm spike, 10 cm spike and 2 cm leaf. Using 
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a ZmGLUlP-GUS construct, Gu et al [68] observed 
a significant up-regulation of the maize /?-glucosidase 
gene, ZmGLUl, in immature transgenic tobacco seeds. 
Vyroubalova et al [33] also detected high expression of 
ZmGLU in immature ears and young leaves of maize. 
Therefore, we speculated that this gene may play 
an important role in reproductive organ initiation and 
seed development. 

Interestingly, and in agreement with the pattern shown 
in our preliminary study [41], at anthesis, when expres- 
sion of IPT, CKX and ZOG genes was uniformly low in 
the flag leaf, the TaGLUl gene showed a marked peak 
of expression, with over a 300-fold up-regulation. From 
anthesis, the flag leaf provides significant resources 
to the developing seed. We postulate that the high p- 
glucosidase gene activity at this stage could lead to the 
release of cytokinin from conjugation in the leaf and 
provide a source of cytokinin for the developing carpels/ 
seeds. Clearly, experiments to support this supposition 
still need to be performed. 

The modest correlation of response regulators with 
TalPT gene expression during the early stage of seed 
development when endogenous cytokinin levels are sub- 
stantially elevated [19-21] indicates that response reg- 
ulators may not be good surrogates of endogenous 
cytokinin content. However, the elevated expression of 
members of all four gene families at 2 daa is indicative 
of dynamic turnover of cytokinins at this key stage of 
seed development. 

An interesting trend shown by our data, and that of 
our preliminary study [41], is that several of the highly 
expressed genes in the seeds have low expression in the 
leaves and vice versa, further suggesting that expression 
of individual family members of the cytokinin synthesis 
and metabolism genes are tissue specific. As a monocar- 
pic plant, wheat development is characterized by a rapid 
translocation of metabolites from leaves to developing 
grains after anthesis. Consequently, delayed leaf senes- 
cence at a late stage of seed development may not be 
beneficial to the final grain yield due to the decreased 
nutrient translocation from leaves to seeds as suggested 
by Sykorova et al [69]. The independent and variable 
expression profiles of different members within or among 
cytokinin regulatory gene families provides useful flexi- 
bility for independently manipulating the endogenous 
bioactive cytokinin levels in reproductive and vegetative 
tissues to achieve maximal seed yield. 

The expression levels shown in this study represent at 
least one of the three homoeologues for each locus in all 
cases apart from the TaZOG2-l, TaZOG2-2, TaGLUl-1 
and TaGLUl-2 genes, in which the three homoeologous 
genes showed different expression profiles. The third 
homoeologue in these cases displayed little activity in 
the initial screen for highly expressed genes. The 



expression among homoeologous genes of other cytoki- 
nin regulatory genes may also vary, or individual homo- 
eologous genes could even be silenced in some or all of 
the tissue samples analysed. Homoeologous gene silen- 
cing and/or differential expression is a phenomenon 
widely documented in bread wheat [43-46] and other 
polyploid species including tetraploid cotton [70], and 
tetraploid Tragopogon mirus [71]. 

Given that homoeologous expression can be differen- 
tially regulated by genetic and epigenetic mechanisms in 
a tissue and developmentally-specific manner [43-45], 
detection and pyramiding of interesting homoeologues 
for improvement of target traits has recently been used 
in wheat breeding programs [72]. For instance, in an 
attempt to improve the bread-making quality of hexa- 
ploid wheat, homoeologous recombination has been used 
to increase the copy number of the stronger flour quality 
contributing homoeologue, Glu-Dl, of the HMW-GS gene. 
A partially isohomoeoallelic line, in which the chromo- 
some 1A homoeologue was replaced by the chromosome 
ID homoeologue, showed significantly improved bread 
making quality. The improved quality can be explained 
by the duplication of the Glu-Dl homoeologue [73,74]. 
Using a different approach, TILLING in combination with 
homoeologue-specific primers, Slade et al [75] obtained 
a triple homozygote loss-of-function of the Wx gene 
and obtained a line with waxy phenotype and improved 
grain quality, 

Therefore, from a breeding perspective, the above 
strategy, of accumulating targeted homoeologues, could 
also be applied to the cytokinin synthesis and metabolic 
gene families. To improve seed yield, higher endogenous 
cytokinin levels could be achieved through tissue- 
specific increased expression of multiple homoeologues 
of the IPT genes TaIPT2, TalPTS and TaIPT8, or 
through decreased homoeologue expression of CKX 
genes TaCKXl and TaCKX2. This could be achieved 
using a number of genetic tools, such as over-expression, 
down-regulation, TILLING and MAS strategies. How- 
ever, while the differential expression profile of each 
homoeologue elucidated in this study could be directly 
used for genetic manipulation towards the improvement 
of the grain yield, precise genome allocation of these 
positive homoeologues is important for further under- 
standing the genetic control, and for facilitating the effi- 
cient accumulation, of these homoeologues within the 
same locus and across different loci. 

Conclusions 

We have demonstrated that the expression patterns of indi- 
vidual members of the TalPTs, TaCKXs, TaZOGs, and 
TaGIUs multigene families were tissue and developmen- 
tally-specific during spike, ovule, seed and flag leaf develop- 
ment. This supports the suggestions that the four multigene 
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families play important roles in maintaining cytokinin 
homeostasis during reproductive development and will con- 
tribute to grain yield and quality. We suggest that artificial 
disturbance of endogenous cytokinin levels at specific stages 
of organ development through targeted manipulation of 
particular family members is a practical possibility, and have 
identified several key candidate genes for genetic improve- 
ment of grain yield in this species. 
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